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A PROBABILITY DENSITY FUNCTION FOR REFLECTED BEAM MOLECULES 


SUMMARY 


A speculative probability density function is constructed on the 
basis of what would be reasonably expected of the probabilities for the 
thermal speed components of molecules that are ejected from an elemental 
surface struck by a molecular beam. It contains five free parameters. 
Three combinations of them occur in the general distribution function 
derivable for molecules moving in the plane of incidence and render 
the law more adaptable than Nocilla's earlier version which provides 
only two. A proper choice of the parameter combinations will usually 
cause it to agree with measured distributions. Pertinent methods are 
developed, and four sample adaptations are described in some detail. 

To determine the five probability parameters themselves the number 
distribution in a second plane must be available and, in addition, the 
heat flux to the surface must be measured. Also one has to stipulate 
that all the incoming molecules and no others are ejecteo from the 
elemental surface. The pressure on it and the shear stress in it can 
then be numerically computed, so that the aerodynamic force acting on 
it is known insofar as it originates with the reflection process. 


I. INTRODUCTION 


A plentiful supply of documentation on the intensity flux measured 
in locally reflected molecule ensembles has long revealed the fact that 
the number of beam molecules returned in the several directions cannot 
always be described, as a simplest approach recommends, by a sphere 
touching the elemental surface imnacted, when it would be proportional 
to the cosine of the angle made by an outgoing direction and the sur- 
face normal. The speeds of such molecules rarely obey a Maxwellian 
equilibrium distribution with perfect accommodation to the surface tem- 
perature (completely diffuse reflection*). Experimental work usually 
operates with a molecular beam generated in a high vacuum chamber and 
often counts only those molecules which leave in the plane of incidence. 


The evidence for largely specular reflection seems more scanty still, 
if net virtually absent. 



If the numbers observed in selected directions are graphically repre- 
sented by proportionate lengths of radius vectors drawn from the point 
of impact (origin), the vector tips occasionally trace out a nearly 
circular curve. In the great majority of cases they point to a roughly 
elliptical ("lobular") shape which is often tilted toward the surface. 
The particular form suggested by a data set depends on the impacting 
material and on that impacted, on the surface structure and contamina- 
tion, on the angle of incidence, and on the energy in the incident 
stream. A rather bewildering variety of forms ensues from these 
concurring conditions. 

A molecular beam will be acting on a surface element when it is 
moving through a rarefied gas. The experimental evidence then indicates 
that the cosine law of reflection is not in general applicable in free 
flight and should not be employed as a matter of course in calculating 
the aerodynamic force attacking a vehicle in the uppermost reaches of 
the atmosphere. 

Since in a beam every molecule possesses, in addition to its random 
speed components, a directed ("macroscopic" or "mass") velocity, it 
seems natural to conjecture that a similar mass velocity may also be 
present after reflection. In reference 1 Nocilla, who introduced that 
idea, was indeed able to approximate, in a satisfactory manner, many 
(but not all) of the experimental point arrays obtained by Hurlbut [2], 
However, he was forced in several instances to postulate the reflected 
mass flow as approaching the surface element, instead of leaving it 
as would seem appropriate. This queer circumstance is apt to throw some 
doubt on the actual existence of a superimposed directed velocity after 
reflection. 

A probability density function for reflected molecules, without 
leaning on that notion and yet without giving up Nocilla' s accomplish- 
ments, can be constructed on the strength of two broad assumptions 
regarding the probabilities for the occurrence of the thermal (random) 
speed components £ 2 » £ 3 . Whether or not the function generated in 

this manner is a solution of Boltzmann's equation is left an open 
question; very likely it is not. But if it can reproduce many experi- 
mental point groupings, especially those not amenable to the interpre- 
tation of reference 1, as in fact it can, one may live with it. It 
may then be considered a bona fide semi- empirical approximation to a 
rigorous solution. 



t 


II. CONSTRUCTION OF THE FUNCTION 


Since Nocilla was able to accommodate a number of experimental 
patterns, it stands to reason that the general form of his approximating 
expression is suitably chosen for obtaining closed curves such as are 
indicated by the observations. It seems therefore advisable to start 
out, as he does, with Maxwell's distribution function excluding, however, 
the mass flow terms. The probability of finding, in any given volume of 
physical space, a thermal velocity in the differential ^inge 


<l t , (i - 1,2,3) 


may then be written as 


Ae 




d£id|2d| 3 , 


(1) 


where A is a constant of proportionality, and 



1 

2RT 


( 2 ) 


(R = specific gas constant; T = gas temperature). With completely dif- 
fuse reflection the molecular velocities are assumed to have adapted to 
the surface ("wall") temperature: T = T w . Expressions (1) and (2) then 

lead to the cosine law which cannot adequately describe the angular dis- 
tribution of re-emitted beam molecules. 


More specifically, the probability of finding a molecule with the 
^-component of the thermal velocity between a fixed value and £ ± + d^ 
emerges as 



x 

1 ? 


-^2 


K I; 


d|. 
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3 




t 


The analogous probabilities for fixed components £ 2 anc ^ £-« arfi 


A — r e 


-* a . 


a 4 d 53 . 


With equal values of £ lt £ 2t £ 3 , these probabilities are equal, a fact 
that characterizes purely diffuse reflection where k 2 is a common con- 
stant. Clearly, one must relinquish this equality in beam reflection 
processes or, phrased differently, adaptation to wall temperature should 
no longer be upheld. Preserving the form of Maxwell's expression as far 
as possible one may define a probability density function as 


AL -(4ti+<fil+*lei> 


(3) 


giving, for fixed values of | lf £ 2 , £ 3 , the probabilities 


A -2- d 5l> A -« 


< 2<3 




d ta . 


K 


1*2 


-< 2 P 2 

e * 3b d| 3 . 


(4) 


which, if the differ, are evidently not equal for equal values of the 
components . 

For a still closer look at the reflection process we first intro- 
duce three physical axes x lf x 2 , x 3 parallel to the axes | lf £ 2 , £ 3 in 
velocity space and, following reference 1, define the beam as arriving 
in the (x 3 ,x 1 )-plane, striking the elemental (x 1 ,x 2 )- - »rface at the 
origin. Any thermal speed component, | 2 , of the reflected molecules 
is normal to that plane and may be expected to be quite as probable as 
its negative counterpart, so that the second of the expressions (4) 
seems adequate. Not the same can be said a priori of the components 
£ x and £ 3 . On the contrary, since some of the incoming molecules 
certainly will not dive into the surface, being repelled by the very 
first surface molecule they encounter, an initial irregularity in the 
distribution of these components arises which is likely to persist even 
if there are subsequent collisions among the reflected molecules." 


In such wall reflection modes, | 2 and -£. still may be considered 
equally likely to occur. 
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Obviously, this irregularity is brought about by something akin to 
specular reflection experienced by a certain fraction of the incoming 
molecules. To account for it and to achieve a difference required by 
it in the probabilities of oppositely equal components, two more prob- 
ability parameters, A x and A 3 , may be employed by writing, in analogy 
to the Maxwellian mass flow terms, (^i-A^" and for and £3 

in the first and third of the expressions (4).* As a consequence, 

£1 * ant * £3 “ A 3 rather than ■ 0 and £3 - Q are the most probable 
values of these components in the reflected gas. 

The modified expression (4) can be set into still another form. 

The probability for a molecule to have any thermal speed is uni^y and 
can most readily be computed from the second formula ( 4 ): 


■boo 


* 3*1 


e-«i*l d* 2 


K. 3 


(5) 


In removing the constant A from the set (4) we may write its modified 
form as 


77 


dfci, 


7! .-«M d is> 

vn 


*2 


d £ 3 . 


( 6 ) 


It is useful to notice that 


> o, 


(7) 


a 3 probabilities are positive quantities. With the method of reference 1 
where 


(2nRT r ) 3 / 2 


1 

v 2RT 

r 


relation (5) is identically satisfied. 
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With surfaces far 
ll into ( 5 2 -X 2 ) 2 . 


from smooth, 


one may be forced alto to alter 


5 



* 


The final expression for the density function proposed becomes 


1 +k,o£2+< 3 ( £3" ^3) J 

n 3/2 


(«) 


The probabilities (6) can be derived from it. With the 

function (8) assumes the form adopted in reference (1) where the quanti- 
ties and X 3 are indeed seen as the components of a reflected mean 
velocity lying in the (x 3 ,x 1 )-plane. The physical oddity occas iona 1 . ly 
encountered with this view did not appear when the angular distribution 
law derivable from the function (8) was applied to four data plots 
selected more or less at random: the mass flow was always directed away 

from the elemental surface, even in two cases where it was approaching it 
in reference 1. Taking this as a general rule, one may at last interpret 
the modifiers and A 3 both ways: either as probability parameters or 

as mass f 1 jw components. 


III. ANGULAR DISTRIBUTION OF REFLECTED MOLECULES 


In deriving the reflection law it is advisable to use spl rical 
cooi Unates in the velocity space with the point of impact at tie crigin 
and the (£ 1 , £ 2 )- plane coinciding with that of the elemental surface, i.e,, 
with the (x 1 ,x 2 )-plane. The azimuthal angle, or longitude, 0 t is counted 
from the positive ^-axis, running from zero to 2n. The co-lacitude, 0, 
counted from the positive | 3 -axla, varies between 0 and n/2 in reflection 
processes, from 0 to n in the reflected gas outside the surface. If v 
denotes the magnitude of the thermal speed, 


£ x = v cos 0 sin 0 
£2 = v sin 0 sin 0 
£ 3 = v cos 0 


where v goes from zero to infinity. The exponential bracket in expres- 
sion (8) assumes the form 


a 2 v 2 - 2bv + c 2 = (av - Q) 2 + c 2 - Q 2 


with 
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(<1 c os ‘‘0 + <§ sin 2 0 ) sln 'e + x 3 cos ! 0 

2 2 

<1^1 c os 0 s In 0 + <3X3 cos 9 

< 1*1 + «§*§ 

b 

a 


(9) 


If dS denotes the area of the elemental surface auu n r is the number 
density in the reflected gas, there will be 

n f dS v cos 0 dt 

molecules of a given speed component S3 returned by the surface during 
the time differential dt. 

The volume element in spuerical coordinates being 


v 2 dv sin 0 d0 60 = v 2 dv dQ, 

the total number of molecules reflected in dt seconds in the direction 
of the solid angle dQ is then indicated by the expression 


n 

r 


K. 

n 3/2 


dS cos 0 dt d(] e 


Q^-c‘ 


e -(av-Q) a V 3 dv 


( 10 ) 


where the probability density function (3) has been employed. The 
evaluation of the improper integral (and others of its kind) can be 
carried out with the aid of reference 3, where a list of basic such 
quadratures is supplied on p. 195/196. One finds its value as 


2 ~T I (l-K) e )e' Q + Q sHT (f + Q 2 ) (1 + erf Q) ] . 


( 11 ) 


The number (10) and therefore the intensity of the molecular flux in the 
elemental solid angle dft is thus proportional to 


C °^4° [(1-K} 2 ) + Q >Tn c| + Q 2 )(l + erf Q)e Q ] 


( 12 ) 
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t 


where 


erf Q 




dx . 


** C 

Note that by the system (9) a depends both on 0 and Q while e is a 

constant which may be taken into the factor of proportionality. 

In order to examine the relationship of expressions (10) and (12) 
with the corresponding results of reference 1, we remark at first that 
there we have 







1 

2RT * 
r 


(13) 


It follows that 


a * K* , b * k 4 (A 1 cos 0 sin e + A 3 cos e) 
b 

Q = “ = k(A x cos 0 sin e + A 3 cos 0). 

a 


Formulas (13) and (16) of reference 1 show that Q is identical with the 
quantity X r which, by the further expressions (14) and (2), is defined as 


X r = s r (sin 0^ cos 0 sin 0 + cos 0 cos 0) 


w j. th 


s 

r 


U 

r 

s 2RT 

r 
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Thus 


A, - U sin e 
x r r 

a 3 ■ u cos e , 
r r * 


(14) 


so that A ! and A 3 indeed emerge as the components of the postulated mass 
velocity U r . Furthermore, 


c = <^Ai + « 3 A 3 - k 2 (A 1 + A§) - < 2 \J 2 * sj. (15) 

-s r 

By this relation and our formula (10) the factor e enters Nocilla's 
expression (13) for the intensity flux. The presence there of the 
denominator nX(o r ) can also be derived from the result (10). To show 
this is somewhat laborious and may be left to the reader if so minded. 

The deviations from the relations of reference (1) center about the 
quantity a which is no longer a constant, but has become a variable 
incorporated, e.g., in expression (12), both explicitly and implicitly 
(by way of Q‘ ). Experimental evidence is as a rule established in the 
plane of incidence where 0=0 (the incoming beam being located in the 
(x 3 ,x i)-plane) . A simplification ensues: 


j a * k\ sin^e + k 3 cos ;? 0 

2 2 

I sin 9 + k 3 A 3 cos e 

V. sIk^ sin 2 0 + < 3 cos -Q 


(16) 


The probability parameter k 2 is absent now in a ; it has no bearing on 
the return in the plane of incidence. The four remaining parameters can 
be reduced to essentially three: 



*3 

\ « 3 Ai(& + l) (17) 

< = *3^3 
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which transform the set (16) into 


< 


a 


<3(1 + 6 cos ’Cp) 


V. 


Q 


ui cos r, + ^ sin r 
v 1 + 6 cos -q? 


( 18 ) 


where at the same time the co-latitude 0 has been replaced by the polar 
angle 


q 



(19) 


in the (x 3 ,x L ;-plane. The second coordinate in this plane, i.e., the 
radius vector, when ascribed a length proportional to the number of 
molecules returned in given directions, will enter into a relationship 
with q, which, on the strength of expressions (12) and (18), may be 
written as 


L = 


(1+& cos ~cp) 


[(1-K) 2 ) + Q 


3 

NMX (“ + 


Q ;? ) (1 + erf Q)e^ ] 


( 20 ) 


where Q is given in the set (18) and the constant I/k^ has been absorbed 
in the factor, 0 , of proportionality. The latter's magnitude is essen- 
tially determined by the scale chosen for plotting the measured flux 
data . 


The function L(q) represents a general angular distribution law 
for molecules reflected in the plane of incidence. It is more flexible 
than the law proposed in reference 1 where the parameter 6 always has 
the value zero. This can be seen by relations (13) and (17). 

Similar but somewhat more complex laws can be written down when 
measurements are to be interpreted, made in planes 0 = 0 C t 0, or, for 
that matter, in planes of constant co-latitude (0 = 0 C ) where the longi- 
tude ;; is variable rather than 0. A plane of the latter kind is con- 
sidered in section V. 


10 
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IV. CURVE ADAPTATION 


To fit the curve L(<p) to a plotted point set, suitable values must 
be found for the parameters 3, 6, " lt ^ of which the last three are 
intimately interlocked in expression (20). As a consequence, the normal 
equations arising with the least square approach are unwieldy to a degree 
rendering the solution for the four parameters all but forbidding except 
perhaps with the use of high-speed computers (which was counter ind icated 
for a number of reasons). Less precise, but better manageable, means were 
substituted for Gauss' method. The adaptations are therefore not the 
best obtainable. It is hoped, notwithstanding, that the four cases 
reported on later do suggest sufficient flexibility of the function L(cp) 
in general, so that the probability density function from which it is 
derived recommends itself as a workable proposition. 

The four specimens dealt with here, and many more of measured point 
arrays, permit us to locate, with some certainty, that particular point 
of the otherwise still unknown curve fit which is farthest away from the 
origin. Quantities related to that maximum of L(<p) are denoted by an 
asterisk. The condition 


±LL^ .iJ 

dcp 


crcp 7 


0 


takes the explicit form 

(1+3 sin 2 cp*) W*‘ + p 3 G* sin cp* • W* - (Q*G*+4)sin 2 cp* = 0 (21) 

where, for any value of cp. 


W = 1 + 6 c os 2 rp 




V. 


G = 


Q (5 + 2Q 2 ) +n/T [| + Q 2 (6 + 2Q 2 ) ] (1 + erf Q)e Q 
(1 + Q 2 ) + Q snr(| + Q 2 ) (1 + erf Q)e Q " 


( 22 ) 
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The function G depends on ; through the variable Q and has the form of a 
logarithmic differential quotient (the numerator being the Q-derivative 
of the denominator '). It grows beyond all bounds when Q -» <», but as can 
be shown with the use of the asymptotic expansion'' of 1 - erf x, tends 
to zero when Q — > -oo. Its values around Q = 0 are moderately large as 
appears from the following table. 


TABLE I 

The Function G in |q| -j 1 


\ 

:q 

-1.0 

-.9 

— 

-.8 

r 

-.7 

— 

-.6 

c 

• • J 

- ■ ■■■■■ 
* 

-.4 

r 

-.3 

-.2 


G 

1.897 

1 .959 

2.024 

2.092 



2.163 

2.237 

2 - 314 

2.395 

2.479 

2.567,2.659; 


Qj 

o ; 

.1 ! .2 

.3 1 .4 

.5 

. 6 

.7 

.8 

.9 

1.0 

G 

2.659 

2 . 754 | 2 .853 

2.955 | 3.061 

3.171 

1 

3.285 

3.48? 

3.523 

3.467 

3.775 


The second differences being nearly equal in this row of G-values, the 
function can be approximated by the quadratic expression 

G « 2.659 + 0.937 + 0. 177 Q ; - 

which in fact is not bad even somewhat beyond [Q = 1. 


One observes that, curiously enough, the denominator itself, which is 
the bracket in expression (20), may be written as 


\ n 

8 

Which advantages 
explored yet. 

** 

See, for example, 


(1 + erf Q)e 
(if any) might be 


Q 


2 


gained 


reference 4, p. 24. 


from this fact has not been 
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Suitable measures of cp* and L* can be read off the given point plot. 
If in addition tentative values are placed both on ^ and Q*, the positive 
solution W* of equation (21) defines the three parameters 


r w* - l 

o — 


cos Cp* 

- . Q*W* - ^ sin q* 

^ 1 cos (p* 


3 - 


, and 


L*W* 4 


(23) 




in cp*[l + Q* 2 + Q* >Tn + Q*)(l + erf Q*)e Q *^] 


7 . 


A table of the error function is supplied in reference 4, p. 24. The 
validity of the two guesses must be checked with further information to 
be drawn from the experimental evidence. After selecting two more 
points (cp,L), preferably far away from the curve's tip on both sides of 
it, one computes the two theoretical values (20) of L from cp and the 
above three parameters ((_i 3 has been chosen beforehand). If they don't 
agree with the experiment. ’ lengths, one has to resort to different 
guesses. By diligently operating in this manner one should be able as 
a rule to reach satisfactory agreement. This trial and error procedure 
has been employed with the last two of our four examples (figures 3 and 4). 
It is surprising and gratifying that, merely using the maximum condition, 
the maximum length, and two further points, one is able to construct a 
rather well adapted curve, a fact that speaks for the suitability of the 
expression L. 

The adaptation in figures 1 and 2 was attained in a somewhat differ- 
ent manner. The given points here permit an estimate of what the rlope 
of the curves should be when they intersect with the vertical axis 
(cp = n/2) . Since 

M-3 = Q Cj) » 

quantities related to that point will be identified by the subscript 3, 
so that, e.g., Q ( tt/ 2 ) = Q 3 . One finds by straightforward differentiation 
that 

tan X, = - 

dcp 

v 2 
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FIG. 1. NITROGEN REFLECTED FROM A CLEAVED 
LITHIUM FLUORIDE CRYSTAL. * - 79 . 9 * 





FIG. 2. NITROGEN REFLECTED FROM A CLEAVED 
LITHIUM FLUORIDE CRYSTAL. 0 L -- -85° 
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where X 3 is the angle of inclination at qp = n/2. Thus a relation between 
parameters is established by the prescribed value of tan X 3 . Only one 
additional point is therefore needed for the checking of guesses; it 
seems natural to decide for the point (n/2, L 3 ) at which the curve ought 
to intersect with the vertical axis. (This choice perhaps should not 
have been made in the case of figure 2 where the left end, although the 
approximation is considerably better than by the model of reference l, 
still is not wholly satisfactory. A point read off nearer to the left 
end might have produced a better result.) 

The importance of equation (21) is evident from the foregoing. Its 
use, however, is precluded when the experimental points do not clearly 
define a maximum distance or when they are arranged symmetrically about 
the x 3 -axis ("upright curves"). In the latter case, ^ = 0* is required 
to ensure equal radius vector lengths L at cp and n - qp, so that 

p .3 s in qp 


Since the maximum distance here is attained at cp = n/2 where W" = 1, the 
value of Q related to it becomes Q* = pi 3 . Equation (21) then reduces to 
the identity 0=0 and is rendered useless. 

In such circumstances it may be feasible to operate with the loca- 
tion (cp, L) of a vertical tangent's point of contact. The general condi- 
tion for infinite slope, 


d dog L ) 

dcp 


'V 

cp-q, 


tan <£, 


leads to the explicit relation 


(1 + 2 sin 2 ^) W + sin q • W - (QG + 4) sin 2 cp = 0 (25) 

which then can substitute for the offices of equation (21). If especially 
the data points group symmetrically about the surface normal, the depend- 
ence of Q on cos cp must be dropped (^ = 0), so that the solution of 


By expressions (17) and (14) this is equivalent to 8 - = 0 in refer- 
ence 1 . 
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equation (25) Is the same whether the angle £* Is used or Its supplement 
n - It should indeed make no difference which one of the two vertical 

tangents one chooses to consider. 

With upright curves there are but three unknown probability param- 
eters (6, ^i 3 , (3). It is tempting to try 6*0 first (Nocllla : s approach). 
With this provision, 


^ 3 H3 sin 


and relation (25) reduces to the equation 


m _ 2 s In £ - 1 
s in cp cos o' 


which can be solved directly for pi 3 because (? depends on (T alone. The 
last remaining parameter, 3 , may be computed, e.g., from the axis length 
L* associated with cp* = n/2: 


0 = 


(1 + 


n!> 


r- ,± 

U-3 V JI (-J 


n!)d 


+ erf (i 3 ) e 


M3 


as Q' = p 3 with a symmetric setup. If the numerical values of and 3 
lead to the proper length L, the assumption 6 = 0 is justified. Other- 
wise, the more cumbersome trial and error procedure involving equation (25) 
must be initiated. 

A spot check of equation (26) was made with the "upright" figure 13 
of reference 1 which is characterized there by 0 r = 0, s r = 1/2. As 
= 0 here, formulas (15) and (17) give 

2 _ 2, 2 _ 2 

s r “ *3^3 “ H3> 


so that (i 3 = 1/2. The above equation is then satisfied with cp = 56.5°. 

The top distance of the figure, L* = 3 cm, gives 3 = 0.701 and L = 1.92 cm. 
Inspection shows that the pair, cp, IT very nearly defines a point of infinite 
slope. As a consequence, curve (13) in reference 1 should be closely 
reproduced by the expression* 


* 

No such representations are given in the referenced article. 
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L ■ 0.701 sin q [1 + Q + Q 


v “ (4 


+ Q ) (1 + erf Q)e ' ], 


where 


-sin 


• • 


All those curves of reference 1 that provide a satisfactory fit 
can he represented in a similar manner.' 1 Our concern is with figures 
(17) ana (18) where Hurlbut's [2] measurements are rather poorly approxi- 
mated. " The results when applying the second method outlined above are 
depicted in our figures l and 2. Major assumptions and the four 
probability part., eters are given in Table II. 


TABLE II 


Significant Quantities for Figures 1 and 2 


1 

, * 

L* 

tan \ 5 

~ — 

3 


hi 

03 

Fig. 1 

2 7° 50' 

7.4 cm 

0.85 

“1 

2.3 cm 

0.398 

-0.605 

0.280 

0.373 

Fig. 2 

27°50' 

7.4 cm 

0.724 

2.5 cm 

0.51 

-0.605 

0.243 

0.324 

i 


The values in the last four columns specify the expressions (20) 
and (18) for L and Q so that the pertaining curves can be plotted. 
Several entries are equal, reflecting a general resemblance of the 
measured point arrays which are both obtained by shooting a molecular 
beam on cleaved lithium fluoride crystals. Merely the beam angles are 
somewhat different (0^ = -79.9° and 0^ = -85°, respectively). The 
almost glancing incidence is perhaps responsible for the strong inclina- 
tion of the data "axis." As broadly indicated by the values of tan X', 
it is even larger in the second case than in the first. 


* 

The not quite commendable adaptations in figures (15) and (16) are 
likely to improve with 5^0. 

One will notice that Nocilla's curve fits generally deteriorate with 
increasing "major axis" inclination which is particularly large in 
these two figures. From this observation it was inferred earlier 
that 5=0 might be a good guess with upright configurations. 
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By the value of 6 we have 



6 + 1 * 0.395 


so that < <-* . Furthermore, expressions (17) and the result (7) show 
that and are both positive. As opposed to reference 1*, the mass 
velocity, if one chooses to interpret the A as their components (14), 
is directed away from the surface, since 0 r has to be an acute angle 
(tan 1.9, or © r ** 62.3° in both figures). 


For further examination of the method two angular distributions were 
selected given in figures 4 and 5 of reference 5 (where the point connec- 
tions shown no doubt have been drawn in with the use of french curves). 
Both depict argon scattering from silver. Since a value of the angle 
X ? cannot be ascertained here, one has to employ the first of the methods 
sketched above. The results are given by our figures 3 and 4. The 
approximations are virtually as good as those by the empirical curves 
excepting two deviating points in the second example which could not 
be reproduced. The abrupt change in curvature they indicate is, however, 
a quite unusual feature attributable perhaps to some freakish surface 
peculiarity. Underlying numerical values are given in Table III. 


TABLE III 


Significant Quantities Related to Figures 3 and 4 



cp* 

L* 

<Pl 

Ll 

Cp 2 



l 2 

3 

6 

Ml 

M3 

Fig. 

3 

50° 

8.3 cm 

CO 

o 

o 

6.0 cm 

o 

O 

5.8 cm 

0.362 

2.13 

2.5597 

0 

Fig. 

4 

o 

00 

8.2 cm 

o 

O 

CM 

4.0 cm 

60° 

3.4 cm 

0.003128 

0 

1.7268 

1.03833 


* 

Here, 6 r is given as 120° in both cases. 

•kit 

Since the angle of incidence is the same (50°), they demonstrate, 
incidentally, that the shape of the reflected curve depends on the 
incident energy (0.064 eV and 4.35 eV, respectively.) Unfortunately, 
free flight energies have not yet been realized experimentally. It 
is most desirable that future efforts correct for this deficiency. 
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The two zero values had suggested themselves in the course of th*» 
tr la 1 -and-error calculations. The first one entails that A, ■ 0 while 


1 ‘ <,(6 + 1 ) 

remains a positive quantity, because » 3 >0 by the result (7). Thus, 
the numerator of the Q-expression contains the cosine term only. Its 
presence causes the shape to be tilted. 

The second zero (A * 0) means that : « «-» and shows that Nocllla's 
approach would do nicely with the distribution of figure 4. Since 
and ^3 are both positive, the reflected mass flow again would turn away 
from the elemental surface. From the formulas (17) one sees that 




tan 0 , 
r 


so that, with the last two figures In the lower row of Table III, the 
angle P r would be 59° equivalent to q. * 31°, while <v* * 38 . Further- 
more, relation (15) shows that, when £ = 0 , 


• r « *3(^1 + A3) = ul + = 4.06. 

Owing perhaps to the slenderness of figure 4, the directed speed of the 
returned molecules thus is surprisingly high; it comes out more than 
twice as large as their most probable thermal speed, whereas in refer- 
ence 1 the ratio was at most unity. 

The above two methods proved to be the relatively most pliable of 
the many one can think of. All of them require four bits of information 
furnished by the raw data. For instance, one may consider the horizontal 
tangent, its point of contact and two more points, or, with tilted curves 
two vertical tangents and their points of contact, or simply four points. 
Whenever possible, it seems advisable to include L* and the maximum condi 
tion ( 21 ). 
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V. THE PROBABILITY PARAMETERS 


The function (20) can be expected to be adaptable to a large number 
of experimental point sets including all that can be handled by the 
method of reference 1. However, the numerical values it yields in a 
given case are too few to completely characterize the probability 
behavior (8) in the reflected gas. 


In the first place, the three parameters (17) that can be elicited 
from data measured in the plane 0=0 evidently do not suffice for the 
determination of K lt , A , A3. A hope to improve on the situation by 
additional observation in a vertical plane 0 = 0 C t 0* or in a horizontal 
plane 0 = 0 C was found abortive. To illustrate this point let us look 
into the second case (where for good resolution we might choose 

0 C = e* = | - <v*). if 


t 



(27) 


the first and last of the expressions (9) may be written as 

r 

a = <3 [ (b + 1) (1 + e sin 2 0) sin :? © c + cos ? 0 c ] 

“ijl cos 0 sin 0c + cos 0 C 

Q - ’ • 

< \(5+l)(l + e. sin 2 0) sin 2 0 c + cos ! 0 C 


(28) 


One gathers from the Q- formula where the parameter combinations and 
P 3 are already known that an approximation of the "horizontal" plot by 
L(.0 will not grant more than a value for (1 4- c ) which like (1 + 6) is 
a ratio, to be sure, this time involving the parameter k 2 for which so 
far no information was available.** 


The symmetry existing about 0=0 with the random speed components 
£ ? and -| 2 taken as equally probable should lead to the same results 
when obser -ing in 0 = -0c» Otherwise, one would have to introduce 
the probability parameter A 2 ; see footnote on p. 5. 


If a satisfactory value of e cannot be determined, the £ -probabilities 
ire not equal on both sides of the plane of incidence and one must 
resort to A 2 ^ 0. The set (9) then becomes more complicated. 


21 



t 


When measurements in a plane 0 = 0 C are evaluated, the expression 
for Q also exhibits but one additional parameter configuration includ- 
ing < 2 * An analogous difficulty evolves with the approach of refer- 
ence 1, in which 6 = e = 0 so that the expression (28) for Q is useless 
for parameter identification and the set (17) reduces to = <A lt 
M3 = <^3 leaving the three quantities Aj, A^ indeterminate. 

To overcome the impasse we will first assume that all the gas 
molecules striking the elemental surface in dt seconds are re-emitted 
and that no others leave the surface. Let us further stipulate that 
the incoming gas obeys Maxwell's velocity distribution law with the 
superimposed macroscopic velocity components 


V L = V t sin e i 
V 3 = V. cos Q^. 


If n^ is the number density, we have the familiar expression 


/~rt 7 _ s 2 

N^tdS = n t J — dtdS [e 3 + s r n; S ,(1 + erf S 3 ) ] 


(29) 


where is the number of molecules striking per unit time and area, 
is the gas temperature, and 


S 3 



With the aid of the probaoility function (8) the corresponding number, 
N r , of molecules reflected per unit time and area becomes 


+00 


-Kn 


N = n 


K 

ei 2 


e -xI(li-Ai) d| 


r 


e“ K ^i d | L , 


^3 e ’" 3d ^3» 


ll =_t 


Ip *- 00 


1 3=0 


or 

n 2 

N — — [e 43 +v~ji 3 (1 + erf n 3 ) ] . 

2 \ r jT <2 


(30) 
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Equating and 


N r yields the ratio 



n^ v 2RT i 


_c- 

e 3 + S 3 (l 4- erf S 3 ) 

2 

t + \ n 43(1 + erf (i 3 ) 


(31) 


which, n r being unknown, still does not specify a value for < 3 . A 
further measurement obviously is called for. That of the heat flux 


q 


e 

r 


(32) 


will be shown to supply the necessary information. 

The internal portions of the energy fluxes follow from the internal 
energy associated with a single molecule which in reference 6, pp. 189/ 
191 is given as* 


1 5 7 Ii 3y raRT 0 3 ) 

where 7 is the ratio of specific heats in the gas. While the incident 
temperature, T^, can be taken as known, the temperature, T r , in the 
reflected gas must be calculated from the mean translational energy in 
it: 


k 3 / 2 


+00 

ff ( 5 f+il+if) e" ^ ^ 1 " ^i> 2+K l6 2+ <l ( £3" ^3) 2 ! d| 1 d| 2 de 3 = 


“ f <2^ + 2^| + 2^ + ^ + 


* 3 

It is written there as — mRTg, where 3 = 


5 - 3y 
3(7-1) * 
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The sum of the first three terms here is equal to the mean value of the 
squared random velocity, 3RT r , so that 




<1 <1 
(i + -4 + -k) 

*2 S 


2 + i + (1 + e) 

6<V 


(34) 


The internal portion of the reflected flux then becomes 


(e r ) 


int. 


= — N 
12 i 


3 - 3 y 2 + 5 + (1 + €V 
7“1 


(35) 


since the assumption N r = was made. Everything is known at right 
excepting the quantity Experimentation in a plane 0 = 0 C is needed 

to obtain a value for t. 


The translational energy flux of the re-emitted molecules follows 
from the integral expression 


(e r ) 


tr. 


m * 1 * 2*3 

= — n 

2 r 


+00 

O P 


2, 2 2, 2 


* 3 / 2 


[*l(£l“ + <2t2 +< 3(t3~ Vi) 


^3=0 


^ 1 1^ £2^ £3. 


After some lengthy evaluation this becomes 


^ e r^ tr . 


m n r _ _ 2 

_ r + Tk^ + A i + A 5) t e + 'f* h3(! + erf h3) ] 


4 Vit K- 


2 <‘. 


2" 

h3 


Since by the formulas (17) 


2 2 

a; + a; = 


_ 2 


"o + 


2 

^3 


1 T " , 2 (6 + l ) 2 ' <| * 
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the expression may be set into the form 


n 

m r 


< e r>tr. = T^T^f p <6.«.5l.W>. "here 


4 V Jt 


VJ 


4+3+(l+Q-~ ^±_ ( +n 2 

2 + &fl (6+1) ^3 


[e + \'^Tp 3 (l+erf p 3 ) ] - 


(36) 


£±L p’^3 
2 


/gain <4 is the only unknown parameter, as n r / k 3 is available through 
relation (31). Thus, we may write 


(e r ) 


tot. 




5 - 3y 2 + ? +(!+£)•*• 
7 -l 3 


+ 


i n 

1 r 


V 7t 



By the quantity (33) and the known expression for the translational 
energy flux of the incident molecules, the total flux of incident energy 
becomes 


< e iW 



where 



n/2RT 1 


is the molecular speed ratio. 

With the use of the measured heat (32) being transferred in unit 
time to the unit area, it follows that 


2 

K l “ 


! I 2+&f( 3 + " rl + +V7„3<1 + erf U3 )] 


"I 


2RT 


i | ^Y 2 + (2+S 2 )e' S 3 + S 3 (| + S 2 )(l 4- erf S 3 ) - qfcnl^RT^ 


-l 


( 37 ) 
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In cancelling Nj the aggregate 


*3 


has been computed here from expressions (31) and (29). The first terms 
in the numerator and denominator vanish for monatomic gases (which have 
no internal degrees of freedom). The quantity <\ has a definite value 
even when q = 0, or is put equal to zero, as one might be tempted to do 
in default of a measurement. 

The numerical values of the remaining four probability parameters 
are found from the definitions (19) and (27): 


<| = <i(i + o 

< l v 1 + 5 

<| = <f(l + 6)- 1 

. ^3 \I 1 + 

A 3 ■ ^ ' Y 

<1 


(38) 


If one measures the molecular fluxes in a plane 0 = 0 C ^ 0 rather 
than in a plane 0 = 6 C in order to arrive at an expression for , some 
quantity t' will have to be defined instead of e and the result (37) 
will assume a different form. 

In concluding we recall that the energy accommodation coefficient 


a 


a 



e 

w 


can be determined from the measured value for q, provided that, when 
operating with polyatomic gases, the wall temperature, T , is known 
which enters through the quantity (33). 
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VI. RE-EMISSION STRESSES 


To calculate the force vector issuing from the action of re-emitted 
molecules on the elemental surface let us first introduce the unit vectors 
_i , j.» k direction of the axes x x , x P , x 3 . The momentum of a single 
molecule, 

m (ih + J.^2 » 

leads to the sum total of the momenta carried away from the unit surface 
in unit time: 


+00 


mn 


K 2.^2^ 3 

<Cl 2 




3“ ^ 3 )* 1 


1 3=0 


d^dlsd^. 

(39) 


Its negative value gives the stresses acting on the surface. 

Evaluation of the i-component yields the shear stress in x 1 -direc tion: 

mn Ai 2 

t = - — “Z — [ e + '/"rt~P3( 1 + erf ^ 3 )!. 

19 2 sf7 < 3 

which, by virtue of relation (31), may be written as 


/ RT^ _ c 2 _ A, 

T i,r = mn tJ T 1 t e 3 + S 3 (1 + erf S 3 )] = - ^ T i (40) 

where tj is the incident stress. 

There is no component in the x 2 -direction (the integral vanishes): 

t =0. (41) 
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This is a consequence of the assumption that | 2 and -£ are equally 
likely to occur which precludes the introduction of a probability 
parameter A- ^ 0 . 

It was pointed out earlier that, with an "upright" data plot, the 
quantity and therefore A x , must be taken as zero. In such an event, 

reflection does not generate a shear stress at all. Nor does incidence 
as a rule, because it should have to be vertical (V x ■ 0 ) for a sym- 
metrical ejection pattern to appear. 

Expression (40) is the total shear stress, t , when A 2 = 0. It 
has been derived with the use of the probability density function ( 8 ) 
and permits the computation of the accommodation coefficient 


a 



(42) 


The numerical value of available from the set (38); by the result 

(7) it is positive with the parameter (which was found so in the 
four cases dealt with here and altogecher might never turn out negative). 
The incident motion (assigned to the quadrant where x 3 > 0 , x x < 0 ) 
defines Oi and V x = sin 0 as negative. The coefficient o will then 
be smaller than unity. 

The ^-component of the integral (39) determines the pressure 


P r 


— + 'Sn (4 + n.3) (1 + erf (j. 3 ) ] 

2 v n S 2 


(43) 


which, being negative, is directed toward the backside of the elemental 
surface. By the ratio (31) and the set (38) 


1 n 
1 r 

2 \~it 


1 + 6 ^ / ^i e ^ + \fl S 3 (l + erf S 3 ) 

<1 1 'I 2 * ^ 3 (l + erf 43 ) 


(43a) 


where < 1 (> 0 ) can be found from expression (37). 
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The pressure accommodation coefficient 



which is sometimes recommended for use can be calculated when the wall 
temperature is known. The representation in terms of speeds, tempera- 
tures, and probability parameters is rather lengthy and is not set down 
here. 


The aerodynamic force exerted by reflected molecules on the sur- 
face element, 


f r = (jLT r + kp r ) dS, 


(44) 


can be resolved into lift and drag components (which refer to a system 
of coordinates usually different from that affixed to the element). In 
using expression (44), it should be kept in mind, however, that outgassing 
or ingassing must be held negligible"’ on account of the condition N r = Nj 
which was a mainstay in deriving major relations, and, secondly, that 
the basic observations were made with reflected molecules hardly subject 
to collisions with incoming molecules and not at all with those reflected 
by adjacent surface elements. The application therefore is restricted 
to free-molecule or nearly f ree-molecule flow until the probability param- 
eters can be obtained in some manner not relying on such observations. 

Even if efforts in this direction should meet with insurmountable diffi- 
culties and must be given up as hopeless, free-molecule reflection, at 
any rate, is better described on the basis of the function (20) and proper 
experimental evidence than by the cosine law. 


VII. CONCLUDING REMARKS 

The tentative probability density function (8) gives rise to the 
angular number distribution law (20) whose analytic form is rather 
involved, so that the trial and error processes proposed for curve 
fitting purposes require lengthy computations often to be repeated 
many times. One may expect to end up with a somewhat more easily tract- 
able expression of the law when the plane 0 = 0 in which most observations 


* 

There is an opinion that in free flight they would be so. 
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arc made is considered the 
system into which the func 
ing through that idea one 
repres enta t ion 


fundamental plane in a cylindrical , 
tion ( 8 ) would have to be transformed, 
indeed arrives at the slightly simpler 


r,y)- 

Carry 


L 


(1 


"I" 4 

+ 6 cos a <p) ’/ 


— Q + (1 + 2Q*’) (1 + erf Q)e Q 
L v n 


(45) 


where and Q again have the meanings (17) and (18). There is a certain 
connection to expression ( 20 ) in that the bracket there is the \”/4 
multiple of the Q-derivative of the bracket here which, incidentally, 
may be written as 


1 d 2 

2 dQ 7 


(1 4- erf Q)e 


i 


The former bracket is essentially the third-order derivative as mentioned 
in the footnote on page 12. It is perhaps worthwhile to inquire more 
deeply into all these relationships; time did not permit doing this here. 

Almost certainly the distribution law (45) will lead to satisfactory 
approximations of measurements performed in the plane of incidence to 
which it is geared. A spot check was made with the data presented in 
figure 13 of reference 1 which, as was explained in section IV, require 
= 0, p! = 0, and ^3 = 1/2. A close fit then ensues with 0 ' = 0.85. 

The different analytic forms notwithstanding, both the laws (20) and 
(45) describe the measured angular distribution adequately with the use 
of the same probability parameters. They ought to do so, since they 
claim to interpret the behavior of the same reflected gas. The formal 
differences are probably brought about by the imprecise handling of 
infinitesimal volumes usually tolerated in dealing with molecular 
ensembles . 

For a complete (and desirable) determination of all the probability 
parameters a survey in a second plane is needed. Moreover, the law (45) 
does not lend itself to an easy comparison with the theoretical results 
in Noc ilia's paper which uses spherical coordinates. For these reasons, 
although set up first, it was abandoned in favor of the law ( 20 ). 
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